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We study the pairing symmetry of the noncentrosymmetric heavy fermion superconductors 
CeRhSi3 and CeIrSi3 under pressures, which are both antiferromagnets at ambient pressure. 
We solve the Eliashberg equation by means of the random phase approximation and find that 
the mixed state of extended s-wave and p-wave rather than the d + f wave state could be 
realized by enhanced antiferromagnetic spin fluctuations. It is elucidated that the gap function 
has line nodes on the Fermi surface and the resulting density of state in the superconducting 
state shows a similar character to that of usual d-wave superconductors, resulting in the NMR 
relaxation rate l/(TiT) that exhibits no coherence peak and behaves like l/(TiT) oc T 2 at low 
temperatures. 
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1. Introduction 

Recent discoveries of heavy fermion superconductors 
without inversion symmetry have attracted much in- 
terest. CePtsSi 1 was first identified and accompanied 
by the subsequent discoveries of CeRhSi3, 2 CelrSia, 4 
Ulr 5 and CeCoGe3. 6 Besides these heavy fermion sys- 
tems, non-heavy fermion materials such as L^PdaB and 
Li2PtsB were also found. 7 fn all these materials, there 
are nonzero potential gradient W averaged in the unit 
cell due to lack of inversion symmetry, which results 
in the anisotropic spin-orbit interaction expressed as 
eh/4m 2 c 2 (k x VV) ■ <r (= aCo ■ tr), where fc is the 
momentum of a particle and ar is Pauli matrices. The 
anisotropic spin-orbit interaction Cq ■ ex, whose general 
form can be determined by a group theoretical argu- 
ment, 8 leads to many interesting phenomena. 9-19 One 
of the outstanding properties is the parity mixing in su- 
perconducting states, 9 ' n ' 13 ' 14 ' 18 ' 19 i.e. the admixture of 
the spin-singlet and triplet states, which are both well de- 
fined in superconductors with inversion symmetry. The 
pairing symmetry in CePtsSi has been studied both the- 
oretically 13, 14, 16, 18, 19 and experimentally 20-24 and it is 
believed that the s + p wave superconducting state is 
realized. Frigeri et al. pointed out that the spin-orbit 
interaction could determine the direction of the d-vector 
as d oc Cq for which the highest transition tempera- 
ture was obtained. A microscopic calculation with the 
detailed structure of the Fermi surface was done 19 and it 
was concluded that s+p wave state is the most probable 
state. 

Among this new kind of compounds, CeRhSi3 and 
CeIrSi3 have many similarities due to the same crys- 
tal structure: qualitatively similar pressure-temperature 
phase diagrams were indeed obtained from resistivity 
measurements. 2,4 They are both in antiferromagnetic 
(AF) ordered states at low pressures, which are driven to 
the superconducting states beyond the critical pressures 
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where the Neel temperature rapidly decreases. The Neel 
temperature at ambient pressure is Tjv = 1.6K (5. OK) for 
CcRhSi3 (CelrSis) and the superconductivity appears in 
wide pressure ranges with approximate maximal tran- 
sition temperature 1.1K at 2.6GPa (1.6K at 2.5GPa), 
respectively. Moreover, the NMR measurements of the 
relaxation rate l/Tf for CeIrSi3 28 suggest the existence 
of the AF spin fluctuations. Also, the neutron scatter- 
ing experiments for CeRhSi3 29 identified the AF ordering 
vectors as Q = (±0.215,0,0.5) and it is concluded that 
the character of the AF is SDW-like. Besides these ex- 
periments, the recent band calculations elucidated that 
the above two compounds have very similar Fermi sur- 
faces 25 ' 26 (FS). These similarities deduced both experi- 
mentally and theoretically motivate us to discuss the su- 
perconductivities of these two compounds in the same 
theoretical framework. 

In this paper, we study the noncentrosymmetric su- 
perconductors CeRhSi3 and CeIrSi3 with particular em- 
phasis on the influence of AF fluctuations to identify the 
pairing symmetry realized in these systems. We also ex- 
amine the properties in the superconducting state, the 
density of states and the NMR relaxation rate which 
characterize the nodal structure of a gap function on the 
FS. 

This paper is organized as follows. In the next sec- 
tion we introduce the model and briefly mention basic 
properties in the non-interacting case. Then in §3, we 
examine possible types of pairing symmetry by means of 
the random phase approximation. In §4, the characteris- 
tic properties in the superconducting state are discussed, 
and a brief summary is given in §5. 

2. Model 

In CeRhSi3 and CeIrSi3, heavy 4/-electrons around 
the Fermi level play important roles for low energy phe- 
nomena, where 4/ electrons come from Ce ions forming 
a body-centered tetragonal(BCT) lattice 2 as shown in 
Fig. 1. For our analysis of the superconductivity, we start 
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Fig. 1. Unit cell of CcRh(Ir)Si3: only Ce sites (open circles) are 
shown for clarity, which form a body centered tetragonal lattice, 
tl, ti and t% are the nearest, the second nearest and the third 
nearest-neighbor hopping integrals, respectively. The ratio of the 
lattice constants is a : b : c = 1 : 1 : 2. 



with the situation that heavy fermions have already been 
formed through hybridizations with conduction electrons 
and are described by an effective Hamiltonian. Although 
these materials have two kinds of the Fermi surfaces 25-27 
apart from the splitting by the spin-orbit interaction, we 
focus on one of them that has a large weight of the total 
density of states. This enables us to simply describe the 
electrons in the materials by the following single band 
model, 



H = 



k i 



(1) 



£ (fc) = (sinfcj,, - sin fc^, 0) , (2) 
£k = — 2t\(cosk x + cos k v ) + 4<2 cos k x cos k y 
—8^3 cos(fc x /2) cos(fc y /2) cos k z — fi, (3) 

where cj^ = (cfcf, Ckif^ are the annihilation (creation) 
operators of the Kramers doublet. The third term of (2) 
is the Rashba-type anisotropic spin-orbit interaction due 
to the lack of inversion symmetry where its coupling 
constant a is estimated to be less than Q.lep (ef is the 
Fermi energy) according to the band calculation. 25 The 
bare Green's function in the normal state is 

ao + (C (k)/\\C (k)\\)-a i 



G°(k) = J2 



-G° T (k), 



T = ± 



G° T (k) 
\Co(k)\\ 



iu>„ - £kr ' 

£ fe + ra||£ (fc)||, 



where k = (iu) n , k). Note that the Green's functions have 
nonzero off-diagonal elements because of the spin-orbit 
interaction. 

With a natural assumption deduced from the above- 
mentioned experimental and theoretical results that the 
AF in these materials are driven by the nesting of 
the Fermi surfaces, we choose the parameters in the 
above Hamiltonian t\,t2,t3 and filling n so that our 




Fig. 



2. Fermi surfaces in our model. The cross sections at k z = 0, 
= 7r/2, k z = 3-7r/4 and k z = n are shown. 
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Fig. 3. Bare susceptibility x°(</) at w = for a = 0. The result 
is obtained for T = 0.04. The left (right) panel is x°(q) on the 
xy plane at k z = n/2 (xz plane at k y = 0) 



model should be consistent with the band calculation 
and the neutron scattering experiment; {ti,t2,t3,n) = 
(1.0,0.475,0.3,1.055). Here, we define t\ as the energy 
unit. For these fixed parameters, the split Fermi sur- 
faces with spin-orbit splitting ~ 2a are shown in Fig. 
2. Our tight-binding model successfully reproduces the 
characteristic features of the Fermi surfaces obtained by 
the first principle band calculations. 25-27 The matrix el- 
ements of the momentum-dependent susceptibility x° (q) 
at U = are expressed as, 
T x ^ 



k)G°(k). 



(4) 



We show the bare susceptibility x°(<z) = X|fi"f(<?) f° r 
several choices of a at u = in Figs. 3 and 4. At 
a = 0, It has peaks around Q 1± <~ (±7r/2, 0, 7r/2), 
Q 2 ± ~ (0) ^ z7r /'^i 7r /2) as shown in Fig. 3. When the spin- 
orbit coupling a is turned on, the anisotropic spin-orbit 
splitting emerges, which suppresses the nesting of the 
FS and thus slightly affects the behavior of \° : the peak 
structure is a little bit smeared, as seen from Fig. 4. The 
peak structure in x° (q) is qualitatively in agreement with 
the neutron scattering experiment, 29 although its profile 
in momentum space is not so sharp in our model. 

3. Pairing symmetry 

In this section, we study the pairing symmetry by solv- 
ing the Eliashberg equation which is expressed as 



AA SlS2 (/c) 



- E 



V,- 



Xk,k') 



k' S 3 S 4 (Ti<T2 
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Fig. 4. Bare susceptibility x°(l) on the X V plane at k z = tt/2 and 
u = for a ^ (T = 0.04). The left (right) panel is for a = 0.1 
(0.2). 



xG° lS3 (fc')G° 2S4 (-A : ')A CTlCT2 (fc')(5) 

where A is the eigenvalue. We use the random phase ap- 
proximation (RPA), which incorporates large AF fluc- 
tuations at the pressures near the AF critical point, to 
evaluate the effective pairing interaction V. In the follow- 
ing discussions, we neglect the normal self energy, which 
may have little influence on determining the pairing sym- 
metry. 30 The effective pairing interaction V evaluated 
within RPA consists of three parts, 



US Sl S3 &S2S4 &81 S2 

rbub 



lad 



L (M')> 

(6) 

where U is the bare Hubbard repulsion and \/ bub 7 y lad 
are calculated by collecting bubble and ladder diagrams, 
respectively. The bubble terms are 

V*f s (k, k') = v h s £ s (k k') v h ±{k + k'), 

Vsgss(k, k ) = v^g^ s (k — k ), 

where 

= -U 2 X ° SSS M/D h u b (q), 

v^sM = u 2 {-x° Sss M - u x Sss Mx s ggs (q) 
+u x ssss{q)x ssss {q))/Dbub{q), 

D huh (q) = (1 + U X ° lm (q)) (1 + U X ° nu (q)) 

-U 2 X° nn (q)x°uu(l) 

and Vjgj S S4 with other spin indices are zero. The ladder 
terms are 

Vssis(k' k ) = u!;sgg(fc — k ) — w'sj g (fc + k ), 
Vsgfs(k' k ) = v £sss(k ~ k ), 
where 

vlf ss (q) = U 2 XsSsS (q)/DiM, 

vlfsM = U 2 ( X °sssM - Ux%ss(Q)X sssM 

+Ux a sgss(q)xssgs(q))/ D i a d{q), 
Aad(g) = (i-f/x? m (5)) (i-u x ° mi (q)) 

-tf 2 X° m (?)X°U T (?) 

and others are zero. Within RPA, only the V^f s terms do 
not conserve the spins of two particles before and after 
scattering, and the parity mixing is driven only by two 



Green's functions which connect the 4-point vertex part 
with the gap function in eq.(5). Thus, the parity mixing 
effect may not be strongly enhanced by U. 

Generally, the gap function A SlS2 (k) is expressed as, 

A (fc) = (A s (iLu n )d a (k)a + A t (iuj n )d(k) ■ er) ier 2 , (7) 

where A s do and A t d are the order parameters for singlet 
and triplet states, which can have non-zero values simul- 
taneously because of the Rashba spin-orbit interaction. 
We denote (d^) _ ^ 3 = (do,d) hereafter. We can deter- 
mine, by solving the above eigenvalue equation (5), the 
symmetry of the gap function and the transition temper- 
ature T c at which the maximum eigenvalue A max reaches 
unity. 

The pairing symmetries of the singlet states {d^ (fe)}r 
for five irreducible representations of C^ v are listed in Ta- 
ble I. Here, we choose {<i r (fc)}r as dF (k) = d^(k)Co(k) 



Irreducible 




representation 


Basis function 

.A,/, N 



^{9xy(x 2 -y 2 )) d,Q 2 (k) = sin 2k x s\n2k y (co$,2k x — cos2fc y ) 
Bl(d x 2_ y 2) 1 (fc) = (cos2fca; — cos2fc J/ ) 

B2 (d xy ) d^ 2 (fc) = sin 2k x sin 2k y 

E (d xz ) a!g (k) = sin k x sin 2k z 



every representation dF(k) = a! r (fc).Co(fc) 



Table I. The irreducible representations of C4u and the basis 
functions. 



for each representation which is considered to be most 
stable in the superconductors with A <C a. 14 ' 16 The 
harmonic wave functions in Table I give the largest con- 
tributions among the functions which belong to a given 
symmetry mainly because of the factor 1/2 in the prop- 
agating vectors Q 12 ±- 

We solve the Eliashberg equation for all symmetries 
and trace each maximum eigenvalue {A r }r with increas- 
ing U at fixed temperature T = 0.04. In the calculation, 
the first Brillouin zone is divided into 16 x 16 x 16 meshes 
and the number of the Matsubara frequencies used is 
512. We have checked that the following results are qual- 
itatively unchanged for 512 Matsubara frequencies and 
(32) 3 fc-meshes. In Fig. 5, X(U, T = 0.04) at a = are 
shown for five irreducible representations. In the case of 
a = 0, the Eliashberg equation (5) is separated into sin- 
glet and triplet parts and solved independently. For sin- 
glet superconductivity, we can see that, among five sym- 
metries, only A^ for d^ 1 — cos2fc z (extended s-wave) 
can reach unity. Other A s ; n are much smaller than A^*, 
and we cannot see any significant difference among them 
in the present calculation. Regarding the triplet part, all 
A tr i are small and none of them can reach unity. 

We remark that, for Ai representation, A^J for 
d Al (k) = dg 1 (k)Co(k) is very small and negative, which 
means that the triplet (p-wave) channel in the effective 
interaction V is weakly repulsive. 

Let us now inspect the effect of anisotropic spin-orbit 
interaction for Ai symmetry, by computing X(U, T = 



4 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



T=0.04 



k z .n/2 



ky-lt/2 



1.2 
1 

0.8 
0.6 
0.4 
0.2 


-0.2 



- A^in — ■ — 




A 2 sin 




B^in * 




B 2 sin a 




Esin 1 







-0.2 



T=0.04 



A,tr1 




A 2 tri — *— 




B,tri »-- 




B 2 tri g 




Etri - -- - - 























Fig. 7. Bare susceptibility x°(<?) = Xf-f-f-[(<?) at a; = and T = 
0.04 with (ti,t 2 ,t 3 ,n,a) = (1.0,0.405,0.3,1.057,0.1). The left 
(right) panel is x°(?) on tne rry plane at fc z = 7r/2 (the 12 plane 
at k y = n/2). 
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Fig. 5. Maximum eigenvalues of the Eliashberg equation for five 
irreducible representations as a function of U with a = and 
T = 0.04. The upper (lower) panel is Aj^ ax for the singlet (triplet) 
gap function. 



Fig. 8. Maximum eigenvalues 
tions with (ti, ti, £3, n, a) = 
0.04. 
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Fig. 6. Maximum eigenvalues A Al for a = 0, 0.1, 0.2 at T = 0.04. 
For a = 0, A is shown only for the singlet gap function, while for 
a 7^ 0, A is for the singlet-triplet mixing state. 

0.04) at a = 0,0.1,0.2. The results are shown in Fig. 
6. Critical values of U, which correspond to the criti- 
cal pressure for the AF transition, increase a little with 
a, reflecting the change of FS that has a tendency to 
suppress x°. Note that, within RPA calculations, the ef- 
fective Coulomb interaction U e g = [7/min(l — Ux°(Q)) 
is a relevant parameter controlling spin fluctuations even 
in the presence of a. We see that A Al decreases with in- 
creasing a mainly due to the suppressed U e s at the same 
value of U and also due to the mixing with the repulsive 
triplet channel. However, it remains larger than unity, 
suggesting the possibility of the Ai-symmetric supercon- 
ductivity. The ratio of the amplitudes of the singlet gap 
function and the triplet gap function for Ai symmetry is 
very small A t /A s < 0.01 for q/0, which means that 
the properties of the superconductivity is characterized 



dominantly by the singlet part. 

Generally, the Ai symmetric gap function has line 
nodes perpendicular to c-axis on the Fermi surface, but 
the averaged value over the FS (d Al )FS has a non-zero 
value. In our model with this gap function, there exist 
line nodes at k z = ±7r/4, ±3tt/4 on the FS. Regarding 
the nodal structure, our Ai gap function d^ 1 is very sim- 
ilar to that for the usual d-wave superconductivity. 

Because the neutron scattering experiment 29 was per- 
formed at ambient pressure, and the true propagat- 
ing vector under pressures is unknown while the super- 
conductivity occurs at high pressures, we also examine 
the pairing symmetry with another set of parameters 
{t 1 ,t 2 ,t 3 ,n) = (1.0,0.405,0.3,1.057) which give maxi- 
mum values of x°(<?) at Q' ~ (0.35-7T, 0.35-7T, 0.57r) as 
shown in Fig. 7. This ordering vector lifts the degen- 
eracy of Q 1 - (0.57r,0,0.57r) and Q 2 - (0, 0.5tt, 0.5tt) 
in the xy components. We show A calculated for five ir- 
reducible representations at a = 0.1 in Fig. 8. In this 
case, we again find that the Ai-symmctric pairing state 
is most likely to appear, which asserts the robustness of 
the stability of the pairing state with the Ai symmetry 
against a slight change of the Fermi surface. Bi sym- 
metry is the second probable candidate because the gap 
function 1 = cos(2fcr) — cos(2fcj,) is favorable with the 
propagating vector Q' x = Q' y ~ 0.35-7T which could lead 
to the sign change d$ 1 (fe) • dp 1 (fe + Q') < on some 
regions {fc G FS}. We calculate {A r } for other sets of 
parameters with which the propagating vector is of the 
form Q — (S, 8, ir/2) or (6, 0, tt/2) and confirmed that 
Ai-symmetric superconductivity is the most probable. 
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Fig. 9. Fermi surface 
(1.0,0.475,0.3,1.05,0). The 




for (ti,t2,t3,n,a) 
gray regions and white 



regions 



correspond to d^ 1 = cos2fc z > 0, d^ 1 < respectively. The 
line nodes of the Ai gap function are located on the boundaries 
between the two regions. The white arrows connecting the 
areas enclosed by contours represent the scattering processes 
associated with the propagating vector Qi^+- 



This is because d 1 = cos 2k z depends only on k z and 
the changes in the x, y components of Q do not affect 
the main scattering processes for the Ai-symmetric su- 
perconductivity. Thus, unless the effects of pressures are 
not restricted to the fluctuations in the xy components 
of Q, it is most stable. 

Let us consider a possible explanation for the stability 
of the superconductivity with this symmetry in the case 
of {t u t2,t 3 ,n,a) = (1.0,0.475,0.3,1.05,0). We think 
that, in our model, taking into account only the scatter- 
ing processes with k ± k = Q t 2 ± gives us intuitive but 
restricted information, because the peak structures of 
Xo (q) are not so sharp and the shape of the FS is compli- 
cated in the 3D momentum space. Nevertheless we try to 
figure out how these scattering processes on the FS con- 
tribute to the realization of the superconductivity with 
the Ai symmetry. Figure 9 shows the FS and the signs of 
the singlet Ai gap function d 1 = cos2fc z ; the gray and 
the white regions correspond to c?q 1 = cos 2k z > and 
dQ 1 < 0, respectively. The FS has a cylinder-like shape 
along z-axis and there exist wide ranges of hot spots. 
Among these regions, the subsets of the FS connected 
to each other via the momentum Q play an important 
role for the superconductivity, when d 1 (fe) • d^ 1 (k ± Q) 
is negative and large. Such spots might be on the sides 
of the cylinder-like FS as shown in Fig. 9 (enclosed by 
contours) and the area of the spots could be large. Thus, 
the scattering processes drawn with white arrows in Fig. 
9 could mediate the superconductivity. 

4. Density of states and NMR relaxation rate 

We now turn to the properties in the superconducting 
state: the density of states p(ui) and the NMR relaxation 
rate 1/T\T with Ai symmetry d^ 1 . They are expressed 



as 



16 



p(uj) 

1 



-_N„(u), 



(8) 



2tt 2Tcosh 2 w 



(|A„H| 2 + |iV a H| 2 )9,) 




Fig. 10. The density of states in the superconducting state with 
A s (0) = 2.0T C , A t (0)/A s (0) = 0.01 and 7 = 0.01A S (0). The 
solid line is for the Ai-symmetric superconducting state and the 
dashed line is for the conventional s + p wave state. 



where 



N n {u) = -J2 lmG °r R ( w + *7, fe) ■ 



kr 



^ImF T 0ii ( W + H,fc). 



kr 



The normal and anomalous Green's functions in the su- 
perconducting state are given by 

iu>„ + £ fe T 



G°r (fe) 



F r °(fc) 



with 



E kT = yJzl + AL, 

A kT = A s (T)d (k)+TA t (T)\\d(k)\\. 

We assume the T dependence of the order parameters 
as A M (T) = A M (0) tanh (l.74y/T c /T - lj , regarding 

A M (0) and the quasiparticle damping factor 7 as fitting 
parameters. In Fig. 10, the density of states for the ex- 
tended s+p wave superconducting state is shown, which 
is compared with that of the conventional s + p wave 
state. In contrast to the conventional s + p-w&vc state, 
the peak of p(u) at u ~ A is largely suppressed and the 
behavior at u) ~ is proportional to w in our system. 
The latter aspect directly follows from the existence of 
line nodes for the Ai-symmetric gap function. As for the 
former effect, one notices that the strong suppression of 
p(uj) makes its profile quite similar to that in usual d- 
wave states. This means that, in CeRhSi3 and CeIrSi3, 
bulk properties which reflect the nodal structure of the 
gap function may show no essential difference from those 
for <i-wave superconductors. 

Next, we proceed to discuss the NMR relaxation rate 
l/(TiT). Since l/(TiT) generally depends on N n (u) and 
N a (uj) not only at ui ~ but also at u) ~ A(0), the exis- 
tence of line nodes does not necessarily mean that 1/TiT 
in the extended s+p wave superconducting state is sim- 
ilar to that in <i-wave states. Figure 11 shows the tem- 
perature dependence of 1/(T\T) normalized by l/(TiT) c 
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conventional s+p — ■ — 




0.1 1 

T/T c 

Fig. 11. NMR relaxation rate 1/TtT as a function of T. The am- 
plitude of the gap function at T = and the quasiparticle damp- 
ing factor are chosen as A s (0) = 2.0T C , A t (0)/A s (0) = 0.01 
and 7 = 0.01A S (0). The line with open circles is for the Ai- 
symmctric superconducting state and the line with black squares 
is for the conventional s + p wave state. 

for the Ai-symmetric superconducting state and the con- 
ventional s+p wave state. By repeating similar calcula- 
tions for several choices of (A(0),7), we find that 1/T\T 
with Ai-symmetric gap function has no coherence peak 
and behaves as 1/(T\T) ~ T 2 , which is characteristic 
of line-node superconductors. This behavior of the NMR 
l/(TiT) is usually typical for dominant d-wave supercon- 
ductivity but not for dominant extended s-wave one. In 
the present system, however, the FS is highly anisotropic 
and the Ai gap function with line nodes could effectively 
behave like a c?-wave gap function on the FS. 

According to the recent NMR experiments for 
CcIrSi3, 28 l/(TiT) exhibits line-node behavior with no 
coherence peak, which seems not contradictory to our 
results. 

5. Summary 

We have studied the pairing symmetry and the na- 
ture of the gap function in the superconducting state 
in the noncentrosymmctric heavy fermion superconduc- 
tors CeRhSi3 and CeIrSi3. Solving the Eliashberg equa- 
tion within RPA, we have found that AF fluctuations 
could mediate the superconductivity with the parity mix- 
ing of the extended s-wave and p-wave states rather 
than the d + f wave state through the Rashba spin- 
orbit interaction. We have confirmed that extended s+p 
wave state is robust against a slight change of the Fermi 
surface under pressure. In the superconducting state, 
the density of states p{i+>) is very similar to that in d- 
wave superconducting states; suppressed p(u> ~ A) and 
p(bj ~ 0) cx bj. Furthermore, the NMR relaxation rate ex- 
hibits \/(TiT) oc T 2 with no coherence peak at T c as in 
the case of the usual d-wave superconductivity. Our re- 
sults suggest a possible understanding of the recent NMR 
experiment within the extended s+p wave state. 

In the present paper, the effects of the normal self- 
energy are not taken into account. According to the re- 
cent experimental observations, the strong-coupling ef- 
fect may be important in CeRhSi3 and CeIrSi3. 31,32 We 
would like to address this issue in the near future. 
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